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Abstract 

We look at the new two-particle correlation function that can provide the 
space-time information about the source of two charged pions. This func- 
tion clearly reflects the chaotic and incoherent nature of the particle emission 
source disturbed by external forces (fields). We argue that such an investiga- 
tion could probe the size and the temperature of the expanded source where 
two pions are emitted in thermal environment. 



1. The Large Hadron Collider (LHC) at CERN is at the stage to provide particle 
physicists with treasure of data. These data allowed a precise measurement of many 
important parameters of modern particle physics in order to test their consistency 
in particular to some theoretical models. The phenomenon of Bose-Einstein correla- 
tions (BEC) is clear and undeniable part of new physics, complicating the quantum 
statistical description of multi-particle final states at finite temperatures. 

Historically, the investigations of particle correlations have been concentrated 
on pion pairs correlations, however have also been applied to more heavy hadrons, 
quarks, protons, and even gauge bosons, photons, Z-bosons etc. (see, e.g., [1-4] and 
the references therein). 

Recently, the first measurements of BEC in proton-proton collisions at y^i = 
0.9 and 2.36 TeV center-of-mass energies recorded by the CMS experiment at the 
LHC have been carried out [5] . In particular, the size of the particle emission region 
r is seen to increase significantly with the particle multiplicity N. However, the 
intercept (or chaoticity) parameter A has the non monotonic behavior with N. The 
dependence of r and A on both iV and yfs, disputed for a long time ( see, e.g., [6]), 
is now very clearly observed in pp collisions [5] . 

This letter gives the extension of BEC phenomenology based of the first principles 
within the Quantum Field Theory at finite temperature, that opens the new horizon 



in experimental study of BEC itself. One of the main parameters is the temperature 
of the particle (emitting) source under the random external forces (fields) influence 
[7-12]- 

By studying BEC of identical particles, it is possible to determine the time scale 
and spatial region over which particles do not have the interactions. Such a surface is 
called as decoupling one. In fact, for an evolving system such as pp collisions, it is not 
really a surface, since at each time there is a spread out surface due to fluctuations in 
the last interactions, and the shape of this surface evolve even in time. The particle 
source is not approximately constant because of energy-momentum conservation 
constraint. 

In this letter, we make an attempt to demonstrate that the problem of properties 
of the genuine interactions can be explored using experimental data collected at the 
LHC. These data can be analyzed through the compared measures of some inclusive 
distributions and final state correlations. 

2. A pair of identical bosons with momenta p\ and P2 and the mass m produced 
incoherently from an extended source will have an enhanced probability C^ipxiP-i) = 
Ni2(pi, P2) /[Ni(pi) ■ N 2 (p2)} to be measured in terms of differential cross section a, 
where 

N 12 {pi,P2) = - , n , n 

to be found close in 4-momentum space when detected simultaneously, as com- 
pared to if they are detected separately with 



In an experiment, one can account the inclusive density ^2(^1,^2) which describes 
the distribution of two particles in Q (the sub- volume of the phase space) irrespective 
of the presence of any other particles 

11 d 2 n 2 

P2{PUP2) 



2\n events dp x dp2 

where ri2 is the number of particles counted in a phase space domain (pi + dpi,p 2 + 
dp2). The multiplicity N normalization stands as 

^2(^1,^2) dp 1 dp 2 = (N(N - S 12 )), 

n 

where (N) is the averaged number of produced particles. Here, 5 12 = for different 
particles, while 812 = 1 in case of identical ones (coming from the same event). 
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On the other hand, the following relation can be used to retrieve the BEC func- 
tion C 2 (Q): 

C (O) - N{Q) (1) 

where N(Q), in general case, is the number of particle pairs in BEC pattern with 

Q = ^-{Px-P2)»-{P1-P2Y = VM 2 -4m 2 . (2) 

In definitions ([!]) and N re f is the number of particle pairs without BEC and 

p m = (u)i,Pi) are four-momenta of produced particles (i = 1, 2); M = J (pi + p 2 )^ 

is the invariant mass of the pair of particles. 

An essential problem is the estimation of the reference distribution N re f(Q) in 
If there are other correlations besides the Bose- Einstein effect, N re f(Q) should 
be replaced by a reference sample corresponding to the two-particle distribution in 
a geometry without BEC. Thus, the expression (HI) represents the ratio between the 
number of particles pairs N(Q) in the real world and the reference sample N re f(Q) 
in the imaginary world. Note, that N re ^ (Q) can not be directly observed in an 
experiment. Different methods are usually applied for the construction of N re f(Q) 
[1] , however all of them have strong restrictions. One of the preferable methods is to 
construct it directly from data, for example, to use the particles pairs from different 
(mixed) events. 

It is commonly assumed that the maximum magnitude of the function C^Q) is 
2 for pi = p 2 if no any distortion and final state interactions are taking into account. 

In general, the shape of C^Q) is model dependent. Its most simple form is 
[13,14]: 

C 2 (Q) = Co ■ [1 + A exp (-Q ■ r) K ] ■ (1 + eQ), (3) 

where Co is the normalization factor, A is the chaoticity strength factor, meaning A = 
1 for fully incoherent and A = for fully coherent sources; the symbol r is often called 
as the "correlation radius", and assumed to be spherical in this parameterization. 
The linear term in ^ accounts for long-range correlations outside the region of 
BEC. For the experimental data fitting it is often used either k — 1 or k — 2. 

Note that the distribution of bosons can be either far from isotropic, usually 
concentrated in some directions, or almost isotropic, and what is important that 
in both cases the particles are under the random chaotic influence or interactions 
caused by other fields in the thermal bath. In the parameterization ^ all of these 
issues should be embedded in A. 
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To enlarge the quantum pattern of particle production process at finite temper- 
ature, we use C 2 (Q) obtained in [2,3] in the following functional form: 

C 2 ( Pl ,p 2 ) ~ £(iV) {1 + Ai(^) e~ A ^ [1 + A 2 (/3) e +A ^ 2 ] } , (4) 

where A p ^ = (pi — p 2 ) M 9? Mt , (Pi ~ V2Y is the smearing smooth dimensionless gen- 
eralized function. dt^ is the nonlocal structure tensor of the space-time size, and 
it defines the domain of emitted particles. £(N) depends on the multiplicity N as 
f (JV) = (iV(iV - 1)}/(N) 2 . The functions Ai(/3) and A 2 (/3) are the measures of the 
strength of BEC between two particles: Ai(/3) = 7(w,/3)/(l + a) 2 , and the correc- 
tion to the coherence function in the brackets of M) is A 2 (/3) = 2a/ ^/^(co, /3). The 
function 7(0;, 0) calls for quantum thermal features of BEC pattern and is defined 
as 

( a\- ( \ n2 (^) 1 \ - ( a\ 1 - U + UJ ' fti\ 

7(",fl = 7(n) = wM n{u , y ^-1^) = ^ W = — > ^ 

where n(u, (3) is the mean value of quantum numbers for Bose-Einstein (statistics) 
particles with the energy uj in the thermal bath with statistical equilibrium at the 
temperature T = 1/(3. The following condition J2f n f(. u > P) = N is evident, where 
the discrete index / stands for the one-particle state /. 

The important function a = a([3,m) entering Ai and A 2 in Q, the measure of 
chaoticity, summarizes our knowledge of other than space-time characteristics of the 
particle emitting source, and it varies from to 00. 

In terms of time-like L , longitudinal L\ and transverse Lj- components of the 
space-time size s/L^, the distribution A p ^ looks like 

A pR -> A pL = (Ap ) 2 Ll + (Ap 1 ) 2 L 2 + (Ap T f L 2 T , (6) 

where L is treated as the measure of the particle emission time, or even it represents 
the interaction strength of outgoing particles. 

Hence, we have introduced a new parameter ^/L%, which defines the region of 
non-vanishing particle density with the space-time extension of the particle emission 
source. Formula Q must be understood in the sense that exp(— A p s){) is a distribu- 
tion that in the limit L — > 00 strictly becomes a 5 - function. For practical using 
with ignoring the energy-momentum dependence of a, one has: 

C 2 (Q) ~ £{N) {l + \M e~ Q2L2 \l + A 2 (/3) e +Q2L ^ 2 } } . (7) 



The parameter L = L((3, m) in formula ([7]) is the measure of the space-time overlap 
between two particles, and the physical meaning of L depends on the fitting of 
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C*2(Q)-function. L can be denned through the evaluation of the root-mean-squared 
momentum Q rms as: 



f °°d\Q\Q* 


'C 2 (Q,P)- 


1 


f?d\Q\ 


~C 2 (Q,/3)-l 





C 2 (Q,(3) 



C 2 (Q,f3) 



where L and Q rm sW) are related to each other by means of 

-i 1/2 



L{f3,m) 



1 + 



1 + 



i 



4a(ftm) 



7p) 
2 



Qrms \P) 



The following restriction \/3/2 < (L ■ Q rm s((3)) < \/3 is evident, where the lower 
bound satisfies to the case a — > (no any distortion in the particle production 
domain), while the upper limit is given by the very strong influence of chaotic 
external fields (forces), a — > oo. The result is rather stable in the wide range of 
variation of a. 

3. We have to emphasize that there are two different scale parameters when the 
BEC are explored [2,3]. One of them is the correlation radius r introduced in ([3]). 
In fact, the latter gives the "pure" size of the particle emission source without the 
influence of the distortion and interaction forces coming from other fields. The other 
(scale) parameter is the scale L = L st of the production particle domain where the 
stochastic, chaotic distortion due to environment is enforced. This stochastic scale 
carries the dependence of the particle mass, the a-coherence degree and what is very 
important - the temperature dependence. 

The question arises: how can BEC be used to determine the effective stochastic 
scale L st and, perhaps, the phase transition? 

At low temperatures (T < n\Jm? + /z 2 , n = l,2,3,...),in the case of two charged 
pions with the mass m, L st has the form (// is the chemical potential): 



L s t(T) 



where 



3a(N) (Ae p ) 2 (m 2 + /x 2 ) 3 / 4 (T/2 vr) 3 / 2 1 + 



(Ae p ) 2 ~ m 2 
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T 



2-^ 

m v m 



+ 



2 
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and Efj, ~ 4.12 MeV been the energy of the muon with the mass in the decay 

It turns out that the scale L st defines the range of stochastic forces. This effect is 
given by a(iV)-coherence degree which can be estimated from the experiment within 
the function (^(Q) as Q close to zero, 6*2(0), at fixed value of mean multiplicity (iV): 

1 + 7 V2( n ) - C 2 {0) + 7 1/4 (n) A /c , 2 (0)[ 7 1 / 2 (n)-2] + 2 

a(N) = — — V — — - . (9) 

C 2 (0) - 1 

The upper limit of 6*2(0) depends on (N) as well as the quantum thermal features 
of BEC pattern given by 7(71): C 2 (0) < £(JV)[1 - 7 1/2 (n)/2]~ 1 . This upper limit is 
restricted by the maximal value of 2 in the ideal case as (JV) — > 00 and 7(71) = 1. 

Note, that for C^Q) - function Q, the limit a — > 00 yields for fully coherent 
sources with small (N), while a — > case stands for fully chaotic (incoherent) 
sources as (N) — > 00. Actually, the increasing of T leads to squeezing of the domain 
of stochastic forces influence, and L st (T = T ) = r at some effective temperature T . 
The higher temperatures, T > T , satisfy to more squeezing effect and at the critical 
temperature T c the scale L st (T = T c ) takes its minimal value. Obviously, T c ~ 
0(200 MeV) defines the phase transition where the chiral symmetry restoration 
will occur. Since in this phase all the masses tend to zero and a — > 0, one should 
expect the sharp expansion of the region with L st (T > T c ) — > 00. 

Looking through the formulas ^ and ^ one can easily find the non-linear in- 
creasing of the linear size of the emitter source with the charged-particle multiplicity. 
This dependence is clearly seen in the CMS experiment (see the bottom panel in 
Fig. 1). 

In terms of the chaoticity parameter A which is related to a(N) through 



a(N) * 1 ~ A + v/I ^ 1 , < A < 1, 
A 

we expect the non-monotonic behavior of A(iV) with N. This pattern has also been 
observed by CMS (see the upper panel in Fig. 1). The following behavior is evident: 
limjv->.oo cn(N) which tends to lim7v-s>oo r (^o) — > 00. 

The qualitative relation between r and L st is the only one we can emphasize in 
order to explain the mass and temperature dependencies of the source size. The 
dependence of the chaoticity a(N) on the minimum scale cut L st and Q rm s can be 
used to define the fit region for different transverse momentum p^- Such a minimum 
cut on L st introduces a lower cutoff on r. 
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Apparatus or analysis effects which 
may result in the sensitivity of external 
random forces influence, may be inves- 
tigated by studying the dependence of 
the correlation functions on the stochas- 
tic scale L st . This effect is expected to 
contribute strongly at small L st (or large 
a and T). 

4. To summarize: the the- 
oretical model for two-particle Bose- 
Einstein correlation function in case of 
two charged pion pairs in pp collisions 
is carried out for the first time. One of 
the main reasons to study BEC at finite 
temperature is the possibility to deter- 
mine the precision with which the source 
size parameter and the strength chaotic- 
ity parameter(s) can be measured at 
particle colliders. 

Such investigations provide an op- 
portunity for probing the temperature 
of the particle production source and the 
details of the external forces chaotic in- 
fluence. 

In this letter, we faced to C^Q, (3) - 
function in which the contribution of N, 
T, a are presented. This differs from the 
methods already used (see, e.g., [1,6]) 
which resemble the traditional way of 

BEC study, however any qualitative interpretations of r, A, e have not been clarified. 

We find that the stochastic scale L st decreases with increasing T slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is 
approached. 

Actually, the experimental measuring of r (in fm) ensures the precise estimation 
of the effective temperature T which is one of the main thermal characters in the 
particle pair emitter source at given a(N) fixed by both C 2 (Q = 0) and (N). T 
is the true temperature in the region of multiparticle production with dimension 
r = L st , because at this temperature it is exactly the creation of two quanta (pions) 



Fig. 1: Values of the A (top) and 
r (bottom) parameters as a function of 
the charged-particle multiplicity in the 
event for combined (dots) and opposite- 
hemisphere, same-charge (open circles) 
reference samples, at 0.9 TeV. The errors 
shown are statistical only. The points are 
placed on the horizontal scale at the aver- 
age of the multiplicity distribution in the 
corresponding bin (borrowed from [5]). 
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occurred in pp collisions, and these particles obey the criterion of BEC. 

Using the CMS data [5], we point out that the effects of random chaotic forces/fields 
and the temperature T are decreasing with y/s, namely, a(N) =1.58 and T = 107 
MeV at v^= 0.9 TeV, while a(N) =1.38 and T = 88 MeV at y/s = 2.36 TeV. 
Acknowledgements: the author would like to thank I. Gorbunov for the help in 
the preparation of the manuscript. 
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